On the practicality of time-optimal two-qubit Hamiltonian simulation 
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What is the time-optimal way of using a set of control Hamiltonians to obtain a desired inter- 
action? Vidal, Hammerer and Cirac [Phys. Rev. Lett. 88 (2002) 237902] have obtained a set 
of powerful results characterizing the time-optimal simulation of a two-qubit quantum gate using 
a fixed interaction Hamiltonian and fast local control over the individual qubits. How practically 
useful are these results? We prove that there are two-qubit Hamiltonians such that time-optimal sim- 
ulation requires infinitely many steps of evolution, each infinitesimally small, and thus is physically 
impractical. A procedure is given to determine which two-qubit Hamiltonians have this property, 
and we show that almost all Hamiltonians do. Finally, we determine some bounds on the penalty 
that must be paid in the simulation time if the number of steps is fixed at a finite number, and show 
that the cost in simulation time is not too great. 

PACS numbers: 03.67.-a,03.67.Lx 
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I. INTRODUCTION 

A central question of quantum information science is to 
determine the minimal time required to perform a quan- 
tum computation using a set of physical resources known 
to be universal for computation. Our understanding of 
what resources are universal for computation is very well- 
developed, and it is known [2^] that when fast local con- 
trol is available, any unitary dynamics capable of generat- 
ing entanglement is universal for computation. However, 
the question of using these resources in a time-optimal 
fashion is, by comparison, understood relatively poorly. 

This paper considers a particular simplified setting, 
that of time-optimal simulation of two-qubit unitaries 
using a fixed interaction Hamiltonian and arbitrary fast 
local control. Arbitrary fast local control means that 
the evolution of the interaction Hamiltonian may be in- 
terrupted by arbitrary single-qubit operations, and that 
these operations take no time to perform. This assump- 
tion corresponds to certain experimental setups where 
single-qubit operations are performed on a much faster 
time scale than joint operations. Hammerer, Vidal and 
Cirac 0, Q have given a construction for this simulation, 
as well as an elegant expression for the minimum achiev- 
able simulation time. 

The simulation scheme of Hammerer et al. uses, in gen- 
eral, an infinite number of steps to achieve time optimal- 
ity. That is, the interaction Hamiltonian is, in general, in- 
terrupted an infinite number of times by local operations, 
and the time between each interruption is infinitesimal. 



A simulation scheme requiring infinitely many time steps 
is not practical for at least two reasons. First, the orig- 
inal premise that local operations can be performed in 
zero time is no longer valid if one must perform infinitely 
many of them. Second, the effects of noise on such a 
simulation will overwhelm the intended coherent dynam- 
ics. The purpose of our paper is to ask, first, whether 
infinitely many time steps are actually required, in gen- 
eral, for time-optimal simulation? We will find that the 
answer is yes. Indeed, we will show that the overwhelm- 
ing majority of two-qubit interaction Hamiltonians have 
this property. Given this, we then address the question 
of determining how close to time optimal a simulation 
can get, given that one demands a simulation using only 
a finite number of time steps. 

The paper is structured as follows. In Sec. [H] we review 
results about two-qubit time-optimal simulation in the 
limit of fast control. In Sec. 11111 we provide a procedure 
for determining which two-qubit Hamiltonians require in- 
finitesimal time steps when used in this setting. Finally, 
in Sec. II VI we quantify the sacrifice that must be made to 
time-optimality when one insists on having a simulation 
using a finite number of time steps. Sec. IVl concludes the 
paper. 
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The purpose of this section is to introduce notation and 
to review some concepts and results associated with time- 
optimal two-qubit simulation in the limit of fast local 
control. We end the section with an introduction to the 
idea of a "lazy" two-qubit Hamiltonian. 
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A. Notation 

Up to rescaling of the ground state energy, an arbitrary 
two-qubit Hamiltonian can be parameterised as follows: 



H 



3 

£ 



(a • a) + (b ■ a) ® I + V Myo-j cr,-, (1) 



where a = (a x ,a y ,a z ) and 6 = (b x ,b y ,b z ) are real 3- 
vectors, M is a 3 by 3 real matrix, and a = (u\, 02, 03) = 
(X, Y", i?) is the vector of Pauli operators. With respect 
to the computational basis {|0), |1)}, the Pauli operators 
are represented by the following matrices: 



X 



1 

1 



Y 



-i 

1 



Z 



1 

-1 



(2) 



When all the entries My are zero we say that the Hamil- 
tonian is local. Otherwise we say that the Hamiltonian 
is nonlocal. We say that a unitary U is local if it can 
be expressed as a tensor product U = A ® B of single- 
qubit unitaries. Otherwise we say the unitary is nonlocal. 
We shall henceforth restrict the single-qubit unitaries to 
be elements of the special unitary group SU(2) (i.e., the 
group of two-by-two unitaries having unit determinant). 



B. Time optimal simulation 

A simulation scheme to approximate an arbitrary two- 
qubit unitary, U, using a fixed Hamiltonian, H, and ar- 
bitrary local unitaries may, without loss of generality, be 
written as follows. 



U 



(A N (g>B N )e 
. . . (At <g> B-Cje 



iHt 



im i(A ®B ), 



(3) 



where the parameters t n are nonnegative. That is, in 
order to achieve the desired dynamics U we can apply H 
as many times as we wish for arbitrary lengths of time, 
interspersed with arbitrary operations on the individual 
qubits. We occasionally refer to Eq. |J3J as being a circuit 
for U. It is worth noting that the assumption that H 
contains no I® I term and that single-qubit unitaries are 
in SU(2), implies that U is in SU(A). These restrictions 
entail no loss in generality, as they simply take advantage 
of the fact that the global phase of a unitary operator is 
irrelevant. 

Corresponding to the simulation Eq. Q is the interac- 
tion time, which we define to be the total time t±+- ■ -+tN 
for which the interaction Hamiltonian is applied. For a 
given U and H , there are many possible circuits each giv- 
ing rise to a simulation of U. Over this range of possible 
circuits for U, there is a corresponding range of values 
for the interaction time. A circuit which achieves the 
minimum interaction time for a given U and H is said 
to be time- optimal. We define Cu(U) to be the mini- 
mum achievable interaction time for simulating U using 



H. Reference [| gives a simple expression for Ch(U), 
in the two-qubit scenario. To discuss this result, we first 
briefly review the canonical form of a two-qubit unitary 
and two-qubit Hamiltonian operator. 



The canonical form of U and H 



For any unitary U € SU (4) there exists a canonical 
decomposition 

U=(d® Dl ) e -^iX«,X+9 2 Y^Y + 6 3 Z^Z) {C2 ^ (4) 

where Ci, D\, C2, and D2 are single-qubit special uni- 
taries, and 9\, 82, and #3 are unique real numbers satis- 
fying 



tt/4 > 0i > 6 2 > |0 3 | > 0. 



(5) 



Although 15 parameters are needed in order to com- 
pletely specify an arbitrary two-qubit unitary U G 
SU(A), the canonical decomposition shows us that the 
nonlocal behaviour of U can be characterised in terms of 
only three parameters, 61, 02 and 63. We call these three 
parameters the canonical-form parameters of U and the 
operator Ug = e - i ^ x ® x + 6 ^ Y+0 ^ z ) the canonical 
form of U. 

The local parts Ci,Di,C2, and D2 of the canonical de- 
composition do not affect the interaction time, as they 
can be trivially included in the first and last steps, 
Aq <g> Bo and An ® Bn, of a simulation. Therefore the 
canonical-form parameters are all we need to know about 
U in order to calculate the minimum required interaction 
time for Eq. J3J). 

How does one calculate the canonical-form parame- 
ters? For completeness, we review the method given in 
Appendix A of 0. 

In the following it will be helpful to take advantage of 
properties of the so-called magic basis 0, 



|1> = 
|2) = 
|3) = 
14) = 



-^=(|01) + |10)) 



V2 

i 

1 
71 



(100) + |H)) 
(|00) -|11)) 
(|01>-|10». 



It is known that, when expressed in the magic basis, 
local two-qubit special unitaries are real, and canonical- 
form unitaries are diagonal. This means that in the magic 
basis the canonical decomposition looks like U = BDS, 
where R and S are real orthogonal matrices, and D is 
diagonal. The diagonal elements of D can be easily writ- 
ten in terms of the canonical-form parameters of U: if 



3 



we define 



951 = 01 + 


6 / 2 - 


f^3 


lp 2 = 01 - 


02 + 


03 


<^3 = -01 + 


02 + 


03 


954 = -0i - 


02- 


03 



then the diagonal elements of D are e * Vl , e lip2 , e tip3 , 
and e- iV4 . Note that Eq. © together with Eq. J5J im- 
plies that 

37T 37T 

— > > <^2 > V>3 > ¥>4 > -T- ^> 

We have that U T U = S T DR T RDS = S T D 2 S, so the 
eigenvalues of U T U are just the squares of the diagonal 
elements of D. That is, 

eig(U T U) = {e- 2iv \e- 2i,fi2 ,e- 2iV3 ,e- 2iVi }. (8) 

To determine the canonical-form parameters for a par- 
ticular U, we first calculate the eigenvalues of U T U 
(where the transpose is taken in the magic basis), then 
derive the cpj via Eq. (JHJ), and finally solve Eqs. A 
word of caution: the task of deriving the tpj from the 
e~ 2%Vi is not as trivial as it may first seem. The problem 
is that, in general, there is no guarantee that the values 
—2ifj will lie in any particular branch of the logarithm 
function. So, naively taking the argument of e~ 2zVj will 
not necessarily give you —2itpj. A relatively simple pro- 
cedure exists to correct for this problem However 
in the context of Sec. IIIII we will see later that taking 
the logarithm of e~ 2%Vi along the standard branch of the 
logarithm function will suffice to evaluate —2ipj. 

Closely related to the canonical form for a two-qubit 
unitary is the canonical form for a two-qubit Hamilto- 
nian. It is discussed in section V.A of 7], where it is 
referred to as the normal form. Given the purely nonlo- 
cal part 

3 

flj = ^ M V a * ® a 3 ( 9 ) 

».j=l 

of a Hamiltonian H, when it is expressed in the form 
Eq. QJ, the canonical form of H is defined to be the 
unique Hermitian operator 

H s = a 1 X®X + a 2 Y® Y + a 3 Z ® Z (10) 

that satisfies 

Hi = {A® B)H s (A t ig) B t ) (11) 

for some local unitary A eg) B, where a% > a 2 > |«3|. 
The existence and uniqueness of this canonical form is 
established in 0, where it is shown that a\, a 2 and |qt3 j 
are the singular values of the matrix M, and sgn(aa) = 
sgn(detM). 

The canonical form of a Hamiltonian H encapsulates 
the nonlocal behaviour of the evolution of H for very 



small time steps. This can be seen as follows. From 
sections III.B and V.A of we can write 

e~ im = (A (8 B)e~ tHst (C ® D) + 0(t 2 ) 

= {A® B )e- lt{aiX ® x+a2Y ® Y+a3Z ® z \C ® D) 
+0(t 2 ) (12) 

for some local unitaries A,B, C, and D. To order t, the 
evolution of H is given by a unitary having canonical- 
form parameters tct\ , ta 2 and ta 3 . 

D. Expression for Ch{U) 

We are almost ready to review the expression for 
Ch{U) given in 0, Before we do so we review the 
concept of special-majorization. Special-majorization de- 
scribes a particular type of partial ordering of three- 
vectors. Its use in 0, Q allows certain results to be de- 
scribed very succinctly. To define special-majorization, it 
is necessary to first introduce the idea of a special- ordered 
three- vector. Given a real vector (3 — (/3i, f3 2 , fa), the 
corresponding special- ordered vector (3 s is defined as fol- 
lows. The absolute value of the components of (3 s are 
given by the absolute value of the components of (3 re- 
arranged in nonincreasing order. That is, |/3|| = 1/3^(^)1 
for the permutation 7r(j) that gives \(3{\ > \(3 2 \ > |/3f|. 
The definition is completed by specifying that (3f and /3| 
are nonnegative, and that /?| has the same sign as the 
product /?i/?2/33- Then, (3 is said to be special-majorized 
by 7 (denoted [3 ~< s 7) if 

P! < it 

Pl+P s 2 -f3l <7i+7 2 S -7 3 S (13) 
(3 S 1 +(3 S 2 +(3 S 3 < 7 f+ 7 | +7 3 s , 

(14) 

where [3 s and 7 s are the special-ordered versions of (3 and 
7- 

Let H be a two-qubit Hamiltonian having canonical 
form Hs and let U be a two-qubit unitary having canon- 
ical form Us- Then Ch(U), the minimum time required 
to simulate U using H , is given by the minimum value of 
ts > such that either 

-< s at s or (15) 

9+(~,0,0)^ s at s (16) 

holds. 



E. The lazy Hamiltonian 

We now introduce the the central concept of a "lazy" 
Hamiltonian. For a given two-qubit Hamiltonian iJ, we 
define a function r{t) as follows: 

r(t)=C H (e-* m ). (17) 
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That is, r(t) is the minimum total time for which the 
Hamiltonian H must be applied, when it is being used 
together with arbitrary local unitaries, to simulate its 
own action e~ zHt . Such a simulation would be of the 
form 

e - lHt = (A N ® B N )e- lHtN (A N -i ® B N - 1 )e- iHt »- 1 
. . . (Ai <g> B^e - "" 1 (A <g> S ). (18) 

The trivial "simulation" having the single step e~ lHt 
has an interaction time of i. Thus the minimum achiev- 
able interaction time will be no greater than t: 

r(t) < t. (19) 

Under what circumstances will r(t) be less than i? It 
turns out that this question is very closely linked to our 
main question: what are the circumstances under which 
a time-optimal simulation will require infinitesimal time 
steps? 

Consider the class of two-qubit Hamiltonians having 
the following property: 

r(t) <t Vt> 0. (20) 

We shall say that a Hamiltonian is "lazy" if it is nonlocal 
and satisfies Eq. (|2*U|) . 

Proposition 1 If a Hamiltonian is lazy, then the time- 
optimal simulation of any two-qubit nonlocal unitary U 
using H requires infinitely many time steps. 

Proof: Suppose there exists a time-optimal simulation 
scheme, of the form of Eq. J3J), for a nonlocal U using a 
lazy Hamiltonian H , where the number of time steps N is 
finite. At least one of the t n must be nonzero, otherwise 
the simulation would be unable to produce nonlocal dy- 
namics. For such a nonzero t n , consider the correspond- 
ing factor e~ lHtn in the simulation. Since H is lazy, there 
exists a simulation for e~ lHt ™ having an interaction time 
less than t n . If we substitute such a simulation for e~ lHtrl 
back into the simulation for U, then the new simulation 
for U now has a lesser interaction time than it did be- 
fore. However this contradicts the assumption that the 
original simulation was time optimal. Hence the premise 
that the original simulation had finite time steps is false. 
Thus we conclude that any lazy Hamiltonian will require 
infinitesimal time steps when used for the time-optimal 
simulation of any nonlocal two-qubit unitary U. I 
To show that a particular H is lazy, it is sufficient to 
show that r(t) < t for all t in some interval (0, e) for any 
positive e. To see this, note that if there is a simulation 
for e~ lHt with interaction time t s < t, then clearly there 
exists a simulation for e ~' lHnt with interaction time t s n, 
for any positive integer n. Thus, r(t) < t implies that 
r(nt) < nt for all positive integers n. So, if r(t) < t for 
all t e (0, e), then r(f) < t for all t > 0. 



III. GENERAL PROCEDURE 

Which two-qubit Hamiltonians are lazy? We have seen 
in the previous section (proposition ^| that lazy two- 
qubit Hamiltonians require infinitely many time steps if 
they are to be used for time-optimal control, and thus 
are impractical. In this section we provide a simple set 
of sufficient conditions for a Hamiltonian to be lazy, ex- 
pressed in terms of the parameters of the Hamiltonian. 
The parameterisation in Eq. is more general than 
it needs to be for this purpose. We can simplify mat- 
ters by using the fact that a Hamiltonian H is lazy if 
and only if (A ® B)H(A ji ® B^) is lazy, where A and 
B are any single-qubit unitaries. This is a consequence 
of the fact that e~ tHt has the same canonical form as 
e -i{A®B)H{A^®B^)t_ Thus ^ without logs of generality we 

choose to only consider Hamiltonians where the purely 
nonlocal part is in canonical form, that is 

3 

H = I ® (a ■ a) + (b ■ a) 1 + ctjQ-j ® aj, (21) 

where a,\ > a-i > {cxsl- 

Recall from subsection III El that we define a Hamil- 
tonian to be lazy if r(t) < t over some interval (0, e). 
Suppose we could find a Taylor series expansion 

r(i)=r(»)+r( 1 >t + T( 2 »( 2 + ... (22) 

for r(t) in the variable t. Thus, 

r(t) - t = t-<°> + (r« - l)t + r^e + ... (23) 

Then, because we can assume t is small, the correspond- 
ing Hamiltonian is lazy if and only if the first nonzero 
item in the list t^°\ (t^ — 1), t^ 2 \ ... is negative. 
Our procedure involves finding expressions for the first 
few items in that list, in terms of the parameters a, b 
and a of the Hamiltonian. We then find the conditions 
under which each expression will be negative. We find 
that in fact = — 1 = = always, and so 
r' 3 ' is the first term that may be negative. Accordingly, 
in the analysis that follows we consider the behaviour of 
r(t) up to order t 3 , so as to arrive at some nontrivial 
conditions for a Hamiltonian being lazy. 

A. Procedure to find the Taylor coefficients of r(t) 

We seek expressions for the Taylor coefficients of r(t), 
namely, t'°), t^ 2 \ and t^ 3 \ The expression for 

Ch(U) involves the canonical- form parameters 9\, 02 and 
#3 of the unitary U. So we first try to find expressions 
for the canonical- form parameters 02 (t) and 93(f) of 

the unitary e~ lHt . From subsection III CI the canonical- 
form parameters can be expressed in terms of parameters 
ipi, . . . , tp4, where 

eig(e- lHTt e- lHt ) = {e~ i2v \ . . . , e"^ 4 }, (24) 
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with the transpose taken in the magic basis. Thus, 

0{t) ={-\ arg(eig(e-^ Tt e-^)) + mr) * (25) 

where the vector of integers n accounts for the ambiguity 
in taking the argument, and where the down-arrow sorts 
in decreasing order so that we are in agreement with the 
ordering of the ipj in Eq. @ ■ However, since we are only 
interested in the behaviour over a small interval t £ [0, e], 
it turns out we can take ft = 0. This can be seen as 
follows. From the discussion at the end of subsection 
III Dl for small t the canonical-form parameters of e~ lHt 
will be small. Thus, the parameters <p\,. . . , ip4 will also 
be small. But this can only be the case when n = 0. (A 
more rigorous proof of this fact is easily deduced from 
the procedure for finding the canonical-form parameters 
described in [|.) Thus, 



0{t) 



arg(eig(e-^ t e~ iHt )) 



(26) 



for t in some interval [0, e]. 
Now, it is possible to write 



„-iH T t -iHt Kit) 

e e = e w , 



(27) 



where K{t) is given by the Campbell-Baker-Hausdorf se- 
ries (for a derivation see, for example, 0) 



1 



K{t) = -it{H 1 + H) + -{-itf^H 1 , H] + 

^t) 3 ([H T ,[H T ,H]] + [H, [H,H T ]}) 

+ ..., (28) 

where [A, B] = AB - BA. Thus, for t e [0, e] we can 
write 

0(t) = (-i ar g(eig(e*«))) (29) 



2i 



eig(K(t)) 



(30) 



From Theorem 1.10 of 9], a normal- valued operator 
function that can be expressed as a power series 



T(t) = T (0) + T {1) t + T (2) t 2 + 



(31) 



has eigenvalues which are holomorphic functions of 
t. Thus, the entries of the vector — ^eig(K(tj) can 
be expressed as holomorphic functions of t, and the 
components of 0{t) are therefore continuous piecewise- 
holomorphic functions of t, over some interval [0, e]. Ap- 
proximating K(t) to some order in t will give the same 
order of approximation for 0{t): 



m + 0(t n )=(--eig(K(t)+0(t n ))) . (32) 



Define K{t) to be the first three terms in the expansion 
of K{t) in Eq. That is, 



K{t) = -it(H J +H) + -(-it 2 )^ 1 , H] 



12 



{-it) i {[H l ,[H 1 ,B\] + [H, [ff,ff J ]]).(33) 



Then, if we define 



we have 



A(i) = --eig(i^)) 



<?(*) = (A(t)) 1 +0^). 



(34) 



(35) 



We can find the first four Taylor coefficients of each com- 
ponent of X(t) in the following way. Each component of 
X(t) satisfies the characteristic equation 

/(*) = det(K(t) + 2i\(t)) = 0. (36) 

We can substitute a Taylor series for X(t): 

f(t) = det(K(t) + 2i(A^ + A^ 1 ^ + ...)) = 0. (37) 

Since Eq. (|37(l must be true for a range of values of t, 
then all coefficients in a Taylor series for f(t) must be 
zero. Finding expressions for the coefficients /w, and 
solving the equations /w = 0, will give us the coeffi- 
cients \(°\ X^\. . .etc. Note that when we wish to find 
a particular term we need only include terms in the 
expansion of A(t) in Eq. (|37|) to order P . 

How can we use the Taylor coefficients of X(t) to find 
the Taylor coefficients of 0{t)1 Does the ordering opera- 
tion in Eq. (|35|l present a difficulty? Not really. Say that 
we knew the Taylor series for each of the four functions 
Ai(t),. . . ,A4(t). We would simply order these vector com- 
ponents with respect to the zero-order Taylor coefficients: 



A£°> > A 2 
>(o) 



(0) 



> Ai 0) > A 



(0) 



If it happened that none of the 



X y -' were equal, that is A^ 0) > A^' ; > A^ > A^ u; , then 
Eq. I|35(l would immediately imply that 



(0) 



,(0) 



AO) 



0{t) = A(t) + 0(i 4 ) 



(38) 



for t in some interval [0,e]. That is, for small t, (pit) is 
equal to X(t) up to order t 3 , where the components of 
X(t) are arranged so the zeroth Taylor coefficients are in 
decreasing order. In the special case where some of the 
zero-order coefficients A^ , A^ , Ag ^ , A^ are equal, then 
we break the tie by considering the first-order coefficients, 
and if those are equal we consider the next highest or- 
der and so-on. In what follows, we will use the ordering 
scheme as described in this paragraph, so that we may 
use Eq. instead of Eq. (J35J). 



6 



B. Procedure to find the Taylor coefficients of the 
components of A(£). 

In this subsection we describe the calculation of the 
Taylor coefficients of the components of X(t) . 

In the magic basis, the Hamiltonian, H, in Eq. (12 1|) 
reads 



H = 

I oi— o 2 +o 3 
ia 3 +ib 3 
ib 2 



ia 2 - 
\ iai+i 



—ia 3 —ib 3 —ia 2 +ib 2 —iai—ibx" 
-01+02+0:3 ia\~ib\ —ia 2 —ib 2 
—iai+ibi — ol\— a 2 — 03 ia 3 —ib 3 



ia 2 +ib 2 



—ia 3 +ib 3 01+02—03, 



(39) 



The components of K (t) in the magic basis are found 
with the aid of the computer algebra system Maple. This 
is straightforward in any of the standard computer alge- 
bra systems, and the specific form is both complex and 
not particularly illuminating, so we will not reproduce 
the components of K (t) here. Note that each component 
is a third-order polynomial in t. 

Next, we find expressions for f^\...,f^\ To find 
f( n >, we evaluate the expression f(t) in Eq. (|37[) using 
Maple, including terms in the series A^ ** + X(l)t + . . . up 
to at least order t n . Then, /( n ) is given by the coefficient 
of the t n term in this expression. Explicit expressions for 
j( n ) will not be given here as they are rather lengthy and 
not illuminating. They are polynomials in the parame- 
ters of the Hamiltonian. 

The next step is to solve the equations f( n > = 0, n = 
0, . . . , 6, via Maple. The results are as follows. 

• Solving /(°) = yields A^ = 0. That is, the zero- 
order Taylor coefficients of each component of A(t) 
are zero. Thus, from Eq. and Eq. © we have 
= 0. 



• Solving /« = 0, = 0, and /( 3 ) 
no new information about the \( n >. 



provides 



• Solving /( 4 ) = yields four solutions, one for each 
component in X(t). We write them in nonincreasing 
order as follows: 



A 



o 1 + a 2 



Ao = Oi — 02 



l 2 
A 4 



-Oi 



A4 = —oti - a 2 



«3 

0:3 
o 2 + o 3 

«3 



(40) 



This gives O^ 1 ' = (01,02,03). 

Solving /< 5 ) = gives A< 2 ) = 0. Thus, 0^ = 0. 

Solving = gives four solutions to A^ 3 ', so 
long as we assume oi > 02 > 03. Each of the 
four solutions for A*- 3 ) correspond to one of the four 
solutions to A^ 1 ) (which were substituted in turn). 



Thus, we are able to correctly associate each of the 
four solutions to a particular component of the or- 
dered vector A(t). For the sake of brevity we will 
not reproduce the expressions for A, . Rather, we 



9 (3) 

a(3) 
7 3 



just provide the resulting expressions for „ 
1 



(3). 



(— oi(a 



b 2 3 ) 



2a 2 a 3 b 3 

(-a 2 (oi + al + b\ + bl) + 2oia 3 6 3 
(-a 3 (oi + a\ + b\ + b\) + 2oia 2 6 2 



2030262) 
2o 3 ai6i) 
2o 2 ai6i) 



(41) 



The special cases oi = o 2 = o 3 , oi = 02 > 03 and 
Oi\ > a 2 = 03 provide different (and rather more compli- 

cated) solutions for the A] compared with above. Ar- 
riving at the solution in these cases requires solving up to 
/( 10 ) = 0. We will not write out these results explicitly. 
Thus we have 



0(t) ^dt + ^t 3 + 0(t 4 ) 



(42) 



for t in some interval [0,e]. 8^ is given in Eq. 141fl . 
except in the special cases noted above. 



Conditions for laziness 



From section II. C of , the expression for Ch (U) takes 
a simpler form when we have 0i + |#3| < f . In this special 
case, Cjj{U) is given by the minimum value of t s such 
that 



(43) 



where again 9 is the vector of canonical-form parameters 
of U and o is the vector of canonical-form parameters of 
H . This special case certainly holds for the canonical- 
form parameters of e~ lHt when t is sufficiently small. In 
this case Eq. (I43f) is equivalent to 



#1 < Otxts 

0i + 2 - 9 3 < (01 + o 2 - 03)^ 
d\ + 2 + 9 3 < (oi + o 2 + o 3 )i s , 

which is equivalent to 

61 

— < t s 
Oil 

+ 8? — 8 3 

< t s 

< U. 



(44) 



Oil 



a 2 
■ 8 2 - 



o 3 

#3 



Q'l + Q'2 



O3 



(45) 



Thus, 



C H {U) 



at\ Oil + Oi 2 — 03 oi + 02 + 03 



(46) 
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Given Eq. l(4*2*)l . we have 

T(t)=t + T^t 3 + 0(t i ) 

for small t, where r^ 3 ) is given by 



(47) 



t^ 3 - 1 = max ■ 



9 (3) 



9 (3) 
7 2 



fl(3) 
7 3 



? (3) 



9 (3) 
7 2 



fl(3) 
7 3 



cti + a 2 + a 3 

(48) 

It is clear that whenever < 0, the Hamiltonian is 
lazy. It is also clear that r' 3 ' is never greater than zero, 
because that would imply r(t) > t, a contradiction. Wc 
find below the solutions (in terms of the parameters of 
the Hamiltonian) for — 0; all Hamiltonians which do 
not belong to this solution set are guaranteed to be lazy. 
Note, however, that the complement of this solution set 
does not entirely characterise the class of lazy Hamilto- 
nians, since there may be Hamiltonians in this set that 
are lazy due to higher-order Taylor coefficients that are 
negative. So our results may not fully characterise the 
set of all lazy Hamiltonians. 
Let 



B x = 








(49) 


B 2 = 








(50) 


B 3 = 


o[ 3 U 


■ *f - 




(51) 



The coefficient is zero if and only if at least one 
of Bj are zero. It is straightforward to show that for 
ct\ > a 2 > as, 



B l 


= 




«2 


= a-3 = h 


= b 3 = 


(52) 


B 2 


= 




a i 


= -bx;a 2 


= -b 2 ;a 3 = b 3 


(53) 


B 3 


= 


44- 


a i 


= bx;a 2 = 


b 2 ;a 3 = b 3 . 


(54) 



We have arrived at the main result of this paper: 

Result 1 Any Hamiltonian of the form of Eq. \21\) for 

which ai > a 2 > a 3 and for which none of the three 
conditions 



1. 



3. 



a 2 = a 3 = b 2 = b 3 = 
a% = -h;a 2 = -b 2 ;a 3 = 
a\ = bi;a 2 = b 2 ; a 3 = b 3 



hold, is lazy. Such Hamiltonians will therefore need to be 
applied infinitely many times when used in a time-optimal 
simulation of a nonlocal two-qubit unitary. 

These conditions obviously make it very easy to gener- 
ate examples of lazy Hamiltonians, and imply that almost 
all two-qubit Hamiltonians are lazy. Note that the spe- 
cial cases a.\ — a 2 > a 3 , a\ > a 2 — a 3 and ct\ — a 2 = a 3 
yield somewhat more complicated conditions for a Hamil- 
tonian to be lazy. These conditions are complex and not 
very illuminating, but can be obtained using techniques 
similar to those described above, so we will not reproduce 
them here. 



IV. USING LAZY HAMILTONIAN IN FINITE 
TIME STEPS 

The results of the previous section show that almost 
all two-qubit Hamiltonians are lazy. This means that, in 
a simulation circuit, infinitesimal time steps must be em- 
ployed to achieve time-optimality. We now show that, 
despite this requirement, if finite time steps are used 
then the corresponding sacrifice of interaction time is not 
very large — only a small relaxation from strict time- 
optimality is required in order to reduce the number of 
time steps to something practical. 

To make our results concrete, we consider the case 
where the unitary being simulated is the controlled-not 
(cnot) gate. Similar conclusions can be reached in the 
general case by following a similar argument to that be- 
low, and making use of the results of ^(J- It can be 
shown that the minimum time for simulating a cnot 
is Cff(CNOT) = where ct\ is the largest canonical- 
form parameter of the interaction Hamiltonian. When 
H is lazy, can we construct a simulation using a finite 
number of time steps such that the total interaction time 
is not much larger than the optimum C^(cnot)? Such 
a scheme is given in whereby an arbitrary nonlocal 
two-qubit unitary U is applied a finite number of times 
together with local unitaries to simulate a cnot. Us- 
ing the scheme in if U has largest canonical- form 
parameter 9\ such that 



n 



IT 



(55) 



is an integer greater than one, then the scheme can be 
used to simulate a CNOT by applying U exactly n times. 
Of course, we are interested in the case when U = e~ lHA , 
that is U is given by the evolution of an interaction 
Hamiltonian over a time A. The total interaction time 
would then be 



t s = nA 



7rA 

4^i(A)- 



From the previous section, the function 
written, for small A, as 



0i (A) = 
Thus, for small A, 

ts = 



aiA + (3) A 3 + O(A 4 



4a i 



6^ 3) A 2 



Ch(cnot) 



0(A 3 ) 
A 2 
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16a( 



0(A 3 



(56) 

l(A) can be 
(57) 

(58) 
(59) 



This shows that to simulate a cnot gate by applying a 
lazy interaction Hamiltonian in a (finite) number of small 
time-steps, then the penalty in the total interaction time, 
as compared with the optimum, is only of order A 2 . 

As an example, consider a specific interaction Hamil- 
tonian H = 0.1X Cg) X + I ® Z '. Using the results of the 
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FIG. 1: Canonical-form parameter 61(A) of the unitary 
e~ iHA , where H = 0.1X ® X + I ® Z. 




FIG. 2: Total interaction time as a function of the number 
of simulation steps, for the simulation of CNOT using the 
Hamiltonian H = Q.1X ® X + I ® Z. 



step (A) is by numerically solving 



0i (A) = 7 1 - 

4n 



(60) 



No solution to Eq. (|60J) exists for n < 8. This can be seen 
from the fact that tt/ (4x7) = 0.112 . . . , which is greater 
than the maximum value that #i(A) takes. For n equal 
to 8 or greater, a corresponding A can be found. Finally, 
the interaction time t s required to simulate the CNOT is 
calculated via Eq. Ij56|l . The results are shown in Fig. 
The dashed line is the optimal time, Cjj (cnot) = 5ir/2. 
The results clearly show a near-optimal simulation with 
relatively small numbers of time steps. For 20 time steps, 
the total interaction time is just 2.8 percent greater than 
the optimal. 

V. CONCLUSIONS 

We have defined a class of "lazy" two-qubit Hamilto- 
nians, those which can simulate themselves faster with 
the aid of fast local control than with uninterrupted evo- 
lution. When a lazy Hamiltonian is used in the time- 
optimal simulation of any nonlocal two-qubit unitary, we 
have shown that the simulation will require an infinite 
number of steps, and thus will be impractical. We have 
derived a simple set of sufficient conditions enabling us to 
prove that a given Hamiltonian is lazy. This set of con- 
ditions implies that almost all two-qubit Hamiltonians 
are lazy. Finally, we have shown that only a rather small 
sacrifice in the simulation time needs to be made in order 
to use a lazy Hamiltonian in a finite-step simulation. 



Acknowledgments 



previous section it can easily be verified that H is lazy. 
The graph of B\ (A) as a function of A is shown in Fig. ^ 
We choose a range of positive integer values of n, and 
for each n we calculate how long the corresponding time 



Thanks to Chris Dawson for helpful discussions and for 
suggesting the name "lazy" Hamiltonian. Thanks also 
to Guifre Vidal for enlightening discussions. HLH and 
MAN enjoyed the hospitality of the Institute for Quan- 
tum Information at the California Institute of Technol- 
ogy, where part of this work was completed. 



[1] K. Hammerer, G. Vidal, and J. I. Cirac, Phys. Rev. A 

66, 062321 (2002), arXiv quant-ph/0205100 
[2] G. Vidal, K. Hammerer, and J. I. Cirac, Phys. Rev. Lett. 

88, 237902 (2002), arXiv |quant-ph/0112168l 
[3] N. Khaneja, R. Brockett, and S. J. Glaser, Phys. Rev. A 

65, 032308 (2001), arXiv |quant -ph/0006114| 
[4] B. Kr aus and J. I. Cirac, P hys. Rev. A 53, 062309 (2001), 

arXiv quant-ph/0011050 
[5] S. Hill and W. K. Wootters, Ph ys. Rev. Lett. 78, 5022 

(1997), arXiv |qiiant -ph/970 304T] 
[6] A. M. Childs, H. L. Haselgrove, and M. A. Nielsen, In 

preparation (2003). 
[7] C. H. Bennett, J. I. Cirac, M. S. Leifer, D. W. Leung, 

N. Linden, S. Popescu, and G. Vidal, Phys. Rev. A 66, 



012305 (2002), arXivj quant-ph /0107035 v2, updated ver- 
sion of [12j. 

[8] R. M. Wilcox, J. Math. Phys. 8, 962 (1967). 

[9] T. Kato, A short introduction to perturbation theory for 

linear operators (Springer- Ver lag, 1982). 
[10] J. Z hang, J. Vafa, S Sastry, and K. B. Whaley, 

arXiv quant-ph/0212109 (2002). 
[11] M. J. Bremner, C. M. Dawson, J. L. Dodd, A. Gilchrist, 

A. W. Harrow, D. Mortimer, M. A. Nielsen, and 

T. J. Osborne, Phys. Rev. Lett. 89, 247902 (2002), 

arXiv quant-ph/0207072 
[12] C. H. Bennett, J. I. Cirac, M. S. Leifer, D. W. 

Leun g, N. Linden, S. Popescu, and G. Vidal, 

arXiv quant-ph/0107035 vl (2001). 



9 



[13] R. K. Brylinski and G. Chen, eds., Mathematics of Quan- 
tum Computation, Computational Mathematics (Chap- 
man & Hall / CRC Press, 2002). 

[14] J. L. Dodd, M. A. Nielsen, M. J. Bremner, and 
R. T. Thew, Phys. Rev. A 65, 040301 (R) (2002), 
arXiv quant-ph/0106064 

[15] P. Wocjan, D! Janzing, and T. Beth, Quan- 
tum Information and Computation 2, 117 (2002), 
arXiv quant-ph/0106077 

[16] W. Diir, G. Vidal, J. I. Cirac, N. Linden, and 
S. P opescu, Phys. Rev . Lett. 87, 137901 (2001), 
arXiv quant-ph/0006034 i'3 and later versions. 

[17] M. A. Nielsen, M. J. Bremner, J. L. Dodd, A. M. Childs, 



and C. M. Dawson, Phy s. Rev. A 66, 022317 (2002), 

arXiv |quant-ph/0109064| 
[18] G. Vidal and J. 1. Cirac, arXiv[quant-ph/0108076 (2001). 
[19] J. L. Brylinski and R. Brylinski, Universal quantum gates 

(2002), chap. II in [l3|, arXiv quant-ph/0108062 
[20] J. A. Jones and E. Knill, Journal of Magnetic Resonance 

141, 322 (1999). 
[21] D. W. Leung, I. L. Chuang, F. Yamaguchi, and Y. Ya- 

mamot o Phys Rev. A 61, 042310 (2000). 
[22] See [lj |m fijj El El El El E3 and references therein; 

see also [2(1 121| for related work. 



